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Abstract. Bose-Einstein correlations in pairs of identical particles were analyzed in e+e^ multihadron 
annihilations at ~91.2 GeV at LEP. The first studies involved identical charged pions and the emitting 
source size was determined. Then the study of charged kaons suggested that the radius depends on the 
mass of the emitted particles. Subsequenty the dependence of the source radius on the event multiplicity 
was analyzed. The study of the correlations in neutral pions and neutral kaons extended these concepts 
to neutral particles. The shape of the source was analyzed in 3 dimensions and was found not to be 
spherically symmetric. In recent studies at LEP the correlations were analyzed in intervals of the average 
pair transverse momentum and of the pair rapidity to study the correlations between the pion production 
points and their momenta (position-momentum correlations). The latest e^e" data are consistent with 
an expanding source. 

1 Introduction 

Bose-Einstein Correlations (BECs) are a quantum machanical phenomenon which manifests in final mul- 
tihadron states as an enhanced probability for identical bosons to be emitted with small relative four 
momentum Q, compared with non identical bosons under similar kinematic conditions [I1I11I3]- From 
the measured effect it is possible to determine the space time dimensions of the boson-emitting source. 
The BEG effect arises from the ambiguity of path between sources and detectors and the requirement to 
symmetrise the wave function of two or more identical bosons. 

In 1954 the radioastronomers R. Hanbury-Brown and R. Q. Twiss proposed a new interferometry 
technique to measure the angular dimension of a star. It required to measure the mixed intensities in 
two radiotelescopes; the dependence of the correlation on the distance between them yielded the angular 
diameter of the astronomical source [T]. G. Goldhaber et al applied the same principle in particle physics, 
in pp annihilations into two identical charged pions, obtaining the radius of the emitting source, ~1 fm 
i- 

The first LEP analyses on BECs concerned identical charged pions, assuming a spherical emission 
source and yielded the size of the source (R~l fm) and the chaoticity parameter Neutral 7r° were 

then considered [5]. Then the study was extended to neutral and charged kaons in order to determine if 
the source radius depends on the mass of the emitted particles [3 [7] . 

Further analyses were performed to establish if the emitting source radius depends on the particle 
multiplicity [8J . Other studies involved the search for BE correlations in multipions . 

BECs were studied in two and three dimensions and one discovered that the emitting source is not 
spherical [TUj. Many studies were made for WW correlations [TT] and also in v interactions [H]. Finally 
they involved the study of expanding sources and trials to determine the emission time |13j . 

We shall make a brief survey of these studies, concentrating finally on the very recent works which 
prove that even in e^e~ collisions one has expanding sources. 
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Figure 1: Scheme of a measurement of BECs. 
a, b are two sources separated by a distance R; 
A, B are two detectors separated by a distance 
L. The emitted particles go from sources to de- 
tectors as a— >A, b— >B or as a— >B, b^A. In as- 
tronomy L<<R, in particle physics L>>R. 



Figure 2: Distribution of two- photon in- 
variant mass, The smooth curves 
are the total Monte Carlo expectation 
(solid line) and the background expecta- 
tion (dashed line). The shaded region is 
the selected window for the tt" signal. 



2 Experimental Procedure 

A detailed description of the OPAL experiment may be found in ref [14 . The most important subdetector 
for BEG studies is the Central Tracking Detector, the Jet Chamber. For tt" studies we needed also the 
barrel electromagnetic calorimeter. A sample of 4.3 million multihadronic events from decays were 
used. A set of quality cuts was applied and one used cuts specific for BEC studies. 

First, the event thrust axis was computed, using tracks with a minimum of 20 hits in the jet chamber, 
a minimum transverse momentum of 150 MeV and a maximum momentum of 65 GeV. Clusters in the 
electromagnetic calorimeter were used for energies exceeding 100 MeV in the barrel or 200 MeV in the 
endcaps. Only events well contained in the detector were accepted, requiring |cos0thrust| < 0.9, where 
^thrust is the polar angle of the thrust axis with respect to the beam axis. Tracks were required to have 
a maximum momentum of 40 GeV and to originate from the interaction vertex. Electron-positron pairs 
from photon conversions were rejected. The selected events contained a minimum of five tracks and were 
reasonably balanced in charge, i.e. — 'n~^\/{n'^y^ + n~^) < 0.4, where njj^ and are the number of 
positive and negative charge tracks, respectively. About 3.7 million events were left after all cuts. All 
charged particle tracks that passed the selections were used, the pion purity being approximately 90%. 

Since in multihadron events more than 90% of the measured tracks are charged pions, the study of 
BEC for like-sign charged pion pairs was usually performed without proper particle identification and 
without purity correction. This choice introduces a small error in the chaoticity parameter A and in the 
radius R of the emitting region. However some analyses were performed with properly identified pions. 
That required some effective cuts on the fraction of the global solid angle acceptance. 

For 7r°7r°, K^K^ and K^K^ correlations, particle identification was necessary [5l[6]. Fig. 2 shows 
the distribution of the two-photon invariant mass and the 7r° event selection. 
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Figure 3: The BEC distribution C'{Q) for 
charged pions vs Qt^^^^- The smooth solid 
curve is the fitted correlation function (the 
excluded regions contain effects from known 
hadron resonances). 



Figure 4: Radius of the emitting region for 
BECs of 2 identical bosons and for FDCs of 2 
identical baryons produced in e+e~ collisions at 
LEP. 



3 BECs from a static source 

BECs in one dimension. The measured BEC function is defined as the ratio C{Q) = p{Q) / Po{Q) , 
were Q is a Lorenz-invariant variable expressed in terms of the two pion four momenta pi and p2 as 
Q^—-{pi - P2)^, p{Q) = {l/N)dN/dQ is the measured Q distribution of the two pions and poiQ) is a 
reference distribution which should contain all the correlations included in p{Q), except BECs. For the 
determination of po{Q), different methods were used: for identical tt+tt^ and tt^tt" one used the tt+tt" 
sample, but also the event mixing reference sample, where pion pairs are formed from pions belonging to 
different events; also a Monte Carlo (MC) reference sample without BECs was used. 

The correlation distribution C{Q) was parametrised using the Fourier transform of the expression for 
a static sphere of emitters with a Gaussian density: 

C{Q) = N{1 + \exp{-R'Q'^)){l + 5Q + eQ^). (1) 

A is the chaoticity parameter, R is the radius of the source, and N a normalization factor. The empirical 
term {1 + SQ + eQ^) accounts for the behaviour of the correlation function at high Q due to any remaining 
long-range correlation. The largest difference among results from different experiments lies in the choice 
of the reference sample: the statistical errors on R is small, but the systematic uncertainty is large. 

Fig 3 shows a typical distribution of C{Q) versus Q; it is relative to a three dimensional analysis, but 
the observed features are typical of all BECs: notice the BEC peak at low Q and the tail at large Q; the 
solid line is a fit to eq. 1, excluding the Q-intervals indicated in the figure, which contain effects from 
known hadron resonances. 

The same analysis was repeated for KK BECs. A similar analysis was performed on Fermi Dirac 
Correlations (FDCs) for identical fermions: in this case there is no peak at small values of Q, but a dip. 
The analysis gives the radius of the emitting regions as shown in Fig. 4. Note that there probably is a 
decrease of R with increasing mass of the emitted identical particles. 

Fig. 5a shows the variation of the emitting radius with the charged multiplicity of the event [5]: there 
is an increase of about 10% of the radius when the multiplicity increases from 10 to 40 charged hadrons 
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Figure 5: a) Increase of the emitting radius 
with increasing event multiplicity and b) de- 
crease of A. 



Figure 6: The distributions in Qtout and Qt^^^^ 
are broader than in Qt- thus Re > Rt^^^^ ~ 



in the final state. This may be related to the number of hadron jets: one has R^jets > Rsjets > R2jets- 
Notice that there is a corresponding decrease of the chaoticity parameter, Fig. 5b. 

In ref. [TU] it was found that there are Stt BECs, that is after removing the effect of 2tt correlations 
on the Stt sample. The present situation is consistent with the relation 

Rstt ~ R27t/V2 (2) 

In ref. [5] it was found that there arc true multiparticlc correlations up to Stt. 
BECs in two and three dimensions. Multidimensional static analyses were performed in 3 dimensions 
using the Longitudinal Center of Mass System (LCMS): the sum of the impulses of the emitted qq pair 
lies in the plane perpendicular to the event axis, defined by the qq direction. The components of the 
3-dimensional distribution in the longitudinal, out and side projections indicate that the last ones are 
larger, see Fig. 6. Thus the longitudinal radius is about 20% larger than the transverse radius: the 
emitting source is ellissoidical, elongated in the qq direction. 

Comparison of BECs in e+e^ and Nucleus-Nucleus collisions. Fig. 7 shows the BEC functions in 
e~^e^ — !■ hadrons and Pb Pb — s- hadrons: note how much narrower is the distribution in Pb Pb collisions: 
the distribution yields a radius i? ~ 6 — 7 fm for the emissions of pion pairs. 

4 Expanding sources 

BECs have been analyzed in Nucleus- Nucleus collisions in order to find evidence for expanding sources 
due to the formation of a quark-gluon deconfined plasma jlS) I16j. Expanding sources may arise in e~^e~ 
collisions because of string fragmentation [T7]. In order to study BECs in non static, expanding sources 
we analyze the correlation functions 
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Figure 7: Comparison of the BE Correlation 
functions for PbPb and e^e~ collisions. 
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Figure 8: Vyk vs pion pair rapidity Y. Ver- 
tical bars include statistical and systematic er- 
rors. YYK=y corresponds to a source expanding 
boost-invariantly. 
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CMC N^^/N^^^ ' 
in bins of the average pair four-momentum with respect to the event thrust direction 



(Pt 



and of the pair rapidity: 



iri = - In 
' ' 2 



{E1+E2) + {PLI +Pl,2) 

{El + E2) - {pis +pia) 



(3) 



(4) 



(5) 



The experimental distributions in dN/d\Y\ and dN/dkt are in good agreement with the distributions 
from the Jetset Monte Carlo. The dependences of C and C on K were studied in three bins of 
(0.0 < |y| < 0.8, 0.8 < |r I < l.e, 1.6 < |y| < 2.4) and five bins of h (0.1 < fct < 0.2 GeV, 0.2 < fct < 0.3, 
0.3 < fct < 0.4, 0.4 < fct < 0.5 and 0.5 < fct < 0.6 GeV). 

Two-dimensional projections of the correlation function C'{Qi, Qt^^^^, Qtout) ^ single bin of \Y\ and 
fct are shown in Fig. [Hk, b. BEG peaks are visible at low Qf, Qt^idc; Qtout ■ 

To extract the spatial and temporal extensions of the pion source from the experimental correlation 
functions, the Bertsch-Pratt (BP) 



C'(Q^Qw„,Qw) = 7V(l-f Ae-('^'«--+^': 
and the Yano-Koonin (YK) 

C'{q,,q,, qo) = 7V(1 + Ae-(«'«?+^'(«'-"*)'«'+'^'(«''~"'')'^o))F(qt, go) 



— w))^^(g,,Q,^^_^^,Q,^^j (6) 



(7) 
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Figure 9; Two-dimensional (a), (b) and one-dimensional (c), (d) and (e) projections of the correlation 
function C"(Qf, Qt.Me'Qtout) for 0.8 < \Y\ < 1.6 and 0.3 <h < 0.4 GeV. Qt,^, < 0.2 GeV in (a), 
Qe < 0.2 GeV in (b). In (c), (d), (e) the projections at low values (<0.2 GeV) of the other variables. 



paramctcrizations were fitted to the measured correlation functions in intervals of kt and In both 
parameterizations, is a normalization factor, A is the degree of incoherence of the pion sources, related 
to the fraction of pairs that interfere. The parameters N and A, whose product determines the size of 
the BEC peak, are significantly (anti)correlated. The two fimctions F{Qi, Qtaidoi Qtaut) ^ (1 + ^longQt + 
etsideQtside + etoutQtout) ^ud F{qt, qe, go) = (1 + StQt + Seqe + Soqo), where a and 6i are free parameters, 
were introduced in Eq. (6) and (7) to take into account residual long-range two-particle correlations due 
to energy and charge conservation. The interpretation of the other parameters in Eq (6), is: 

- iitside ^^'^ -Riong are the transverse and longitudinal radii in the longitudinal rest frame of the pair; 

- i?tout and the cross-term iiiong.tout are a combination of both spatial and temporal extentions of the 
source. The difference (^tom ~ ^Ldo) proportional to the duration of the particle emission process, 
and i?iong,tout to the source velocity with respect to the pair rest frame. 

In the YK frame Eq. (7), where 7 = — v'^, the free parameters are interpreted as follows: 

- V is the longitudinal velocity, in units of c, of the source element in the CMS frame; 

- Rq measures the time interval, times c, during which particles are emitted, in the rest frame of the 
emitter (source element). The limited phase-space available limits the analysis for Rq ; 

- Rt and Re are the transverse and longitudinal radii, in the rest frame of the emitter. 
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The parameters Rq, Rt and Rg are evaluated in the rest frame of the source element. The two 
parameterizations are not independent, so that a comparison between the BP and YK fits is an important 
test. In the YK picture the source velocity v of each element does not depend on kt, while it is correlated 
with the pair rapidity Y. Yyk measures the rapidity of the source element with respect to the cm 
frame Yyk = 5 In [(1 + — v)]. A non expanding source corresponds to Yyk — for any Y. For a 
longitudinally boost invariant source (for which the velocity of each element is = z/t, where t is the time 
elapsed since the collision and z is the longitudinal coordinate of the element) the correlation Yyk = Y 
is expected as in Fig. 8. 

The following relations hold between the BP and YK parameters: 

Rl.. = Rt (8) 

-Rlong = 7lCMs(^£ + /^LCMS-^o) (9) 

{Rl^^-RlJ = [3hlcMs{Rl + f3lcMsRi)- (10) 

/3lcms is the velocity of the source element in the LCMS, i.e. with respect to the pair longitudinal rest 
frame; 7lcms ~ l/-\/l — PIqms- ^'^^ ^ boost-invariant source, /3lcms — 0: (9) and (10) become: 

-Riong - (11) 
{Rl^^-RlJ PlRl (12) 

In Fig. 10 the best-fit BP and YK parameters are compared: 

- The longitudinal parameter Rf^^^ is larger than i?| in all rapidity intervals. Fig. 9(a),(d) and (g). 
R\ong > Rj corresponds to /Slcms >0, i- e. to a pion source whose expansion is not exactly boost- 
invariant. 

- The equality of the transverse parameters, ^^t^idc i is confirmed; there may be deviations at low kt. 

- Rq and (i?tout ^ Rt^idJ essentially equal to zero, suggesting that the present technique does not 
allow to mesure the duration of the emission process. 



5 Conclusions 

We have first summarized the results obtained on BECs in e+e^ collisions at the Z'^ peak assuming a 
static source. Then we presented an analysis in bins of the average 4-momentum of the pair. Based on 
this, the dynamic features of the pion emitting source were investigated in the YK and BP formalisms. 

The transverse and longitudinal radii of the pion sources decrease for increasing kt, indicating the 
presence of correlations between the particle production points and their momenta. The YK rapidity 
scales with the pair rapidity, in agreement with a nearly boost-invariant expansion of the pion source. 
Phase space limitation did not allow the measurement of the duration of the particle emission process. 

Similar results have been observed in more complex systems, such as the pion sources created in pp 
and heavy- ion collisions, which are now complemented with measurements in the simpler hadronic system 
formed in e+e^ annihilations. The unexplained similarities between BECs in different reactions might 
indicate a present limitation of our understanding of these correlations |18j . 
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